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Abstract 

We use Boulware's Hamiltonian formalism of quadratic gravity theo- 
ries in order to analyze the classical behaviour of Bianchi cosmolog- 
ical models for a Lagrangian density JC = R + f3cB? in four space- 
time dimensions. For this purpose we define a canonical transfor- 
mation which leads to a clear distinction between two main variants 
of the general quadratic theory, i.e. ioY C = R + (3cR? or conformal 
C = oicC"^^'^Caf3iJLv Lagrangian densities. In this paper we restrict the 
study to the first variant. For Bianchi-typc I and IX models, we give 
the explicit forms of the super-Hamiltonian constraint, of the ADM 
Hamiltonian density and of the corresponding canonical equations. In 
the case of a pure quadratic theory C = PcR^, we solve them ana- 
lytically for Bianchi I model. For Bianchi-type IX model, we reduce 
the first-order equations of the Hamiltonian system to three coupled 
second-order equations for the true physical degrees of freedom. This 
discussion is extended to isotropic FLRW models. 

PACS numbers: 0420F, 0420J, 0450, 9880H 

Running title: Bianchi models in quadratic theories of gravity 
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1 Introduction 



Recent propositions of unified field theories have stressed the importance of 
including in the basic Lagrangian, in addition to the Einstein-Hilbert scalar 
curvature term of General Relativity, non linear terms in the various cur- 
vature tensors, leading to the so-called Higher Derivative Gravity Theories. 
The basic motivation for studying these theories comes from the fact that 
they provide one possible approach to an as yet unknown quantum theory of 
gravity. Moreover these theories occur as a necessary low-energy outcome of 
superstring theory even if the issue of the specific choice of the non linear 
terms that have to be kept into the Lagrangian remains somewhat controver- 
sial . Be that as it may, in the very early universe the structure of classical 
solutions of higher derivative gravity may provide a better approximation to 
some low-energy limit of quantum cosmological solutions than those provided 
by General Relativity. Therefore it is of interest to study also the classical so- 
lution space of these theories. In particular, quadratic theory which has been 
studied in great detail is generally viewed as a possible basis for quantum 
gravity f^, since it is renormalizable and asymptotically free |@]. In four di- 
mensions, due to the Gauss-Bonnet topological invariant, the corresponding 
Lagrangian contains only two quadratic terms among three possible ones : 
R^, RapR^^ and Rap^uR^^^^ i where i?, R^p and Ra/Sfiu denote as usual the 
scalar curvature, the Ricci and Riemann tensors respectively. Therefore, the 
i?-squared term as well as a quadratic term in the Weyl tensor, Cap^iuC"''^'^'^ , 
can be chosen as the two independent quadratic terms of the effective action 
of the theory. 

Cosmological applications of this formalism have been considered mainly 
in the context of inflationary cosmology, since the i?-squared term has the 
virtue of inducing an early inflationary stage in the spatially homogeneous 
and isotropic Friedmann-Lematre- Robertson- Walker (FLRW) models [^]. At 
the quantum level. Hawking and Luttrell |0 have examined the effect of 



3 



curvature- squared terms on the wave function of the closed FLRW model, 
in the framework of Hartle-Hawking's path integral approach to quantum 
gravity. However, due to the conformal flatness of FLRW metric, the Weyl- 
squared term plays no role, while the i?-squared one behaves like a mas- 
sive scalar field. For pure fourth-order {C ~ B?) quantum cosmology, the 
Wheeler-De Witt equation has been solved approximatively by Kasper [0 
and exactly for the same model by Pimentel and Obregon . The quantum 
cosmological consequences of introducing a cubic term in the scalar curvature 
into the Einstein-Hilbert action have been investigated by van Elst, Lidsey 
and Tavakol p. 

Although some work has been devoted to the study of the classical dynam- 
ics of Bianchi cosmological models in quadratic theories of gravity, it has 
been performed without any explicit reference to an Hamiltonian formula- 
tion. However any canonical method of quantization of constrained systems 
is based, at a fundamental level, on Hamiltonian concepts [1^ and especially 
on the explicit use of the so-called super- Hamiltonian and supermomenta 
constraints. 

Hamiltonian formulation of General Relativity has been mainly developed in 
classical work by Dirac |[Tl[] and Arnowitt, Deser and Misner (ADM) [|12|. It 
was applied to classical cosmological models including Bianchi ones, giving 
rise to the so-called Hamiltonian Cosmology [l^. Its quantum counterpart. 



Quantum Cosmology [1^1 , is based on the application of canonical or path 
integral quantization methods of General Relativity to the very early stages 
of the universe, i.e. in the neighbourhood of Planck's time. 
Due to the complexity of the task, Hamiltonian formulations of theories with 
non linear Lagrangians in the curvature have up to now been restricted 
to quadratic Lagrangians only. They have been developed first of all for 



four dimensional space-times by Boulware [|1^] and Buchbinder-Lyakhovich 
(B & L) []TB[, and for multidimensional space-times by the two last authors 



and Karataeva '^Tj. The main objective of these formulations was the build- 
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ing of a consistent quantum theory for higher derivative theories of gravity. 
No systematic apphcation of this formahsm to cosmological models at the 
classical or quantum level has been attempted, with the exception already 
mentioned above concerning FLRW models. 

The plan of the paper is as follows. In section 2 we summarize the main 
results of the Hamiltonian formulation of general quadratic gravity theories 
and settle down our notation. In section 3 we use this Hamiltonian formal- 
ism in the case of the Einstein- Hilbert Lagrangian containing a quadratic 
correction term for the classical study of spatially homogeneous cosmological 
models. We give the explicit form of the super-Hamiltonian constraint, of 
the ADM Hamiltonian density and of the corresponding canonical equations 
for Bianchi-type 1 and IX models. In the particular case of a pure quadratic 
theory (£ = PcR"^), we solve them analytically for Bianchi I model. For 
Bianchi-type IX model, we reduce the first-order equations of the Hamil- 
tonian system to three coupled second-order equations for the true physical 
degrees of freedom. In section 4 this discussion is extended to isotropic FLRW 
cosmological models. In the last section we discuss our results and further 
investigations initiated by this work. 
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2 Hamiltonian formalism 



In classical mechanics, it is possible to construct a well defined Hamiltonian 
formalism only if the Lagrangian is regular (i.e. the Hessian determinant of 
the Poisson-Hamilton transformation is not equal to zero). In the opposite 
case, the system is said constrained and we have to use Dirac's approach [jlO 



in order to define a consistent Legendre transformation. When dealing with 
higher order field theories (in which Lagrangians depend on derivatives of the 
generalized coordinates of any order), it is necessary to follow a generalized 



Ostrogradski's method to perform such a transformation [16|. The result- 
ing canonical theory is consistent and it has been proved by generalization 



of Ostrogradski's theorem |]T8l that the Euler-Lagrange equations (derived 
from the Lagrangian by the usual variational principle) and Dirac-Hamilton 
canonical equations (obtained from the explicit form of the total Hamilto- 
nian) are in fact equivalent. 

As mentioned in the first section, Hamiltonian formulation of gravity theo- 
ries with quadratic Lagrangians was achieved by Boulware and B & L, for 
four- dimensional space-times. In this section we give a short summary of 
their main results and in particular we write down the general expressions 
for the super-Hamiltonian, the supermomenta and the canonical equations 
which are our starting point for cosmological applications. 



First of all, we settle down the notations and conventions used throughout 
this paper. We take as the general form of the gravitational action in four 
dimensions the following expression written in terms of scalar curvature and 
Weyl tensorQ : 

where ac and (3c are coupling constants, A is the cosmological constant, G is 
Newton's constant and Cap^iv is the Weyl tensor in four dimensions defined 

^We use MTW conventions 
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by 



I3fiu 



+ R{5lgp,-5:gp,)/Q 



As in ADM method, we assume a foliation of space-time in spacelike 3- 
dimensional hypersurfaces and we use ADM coordinate basis to compute 
the components of the various curvature tensors. The "0" index refers to 
the normal component to the hypersurfaces. Subscripts '•^^ and indicate 
quantities respectively defined over the 4-dimensional space-time and over 
the space-like hypersurfaces. Most of the time, when there is no ambiguity, 
the last subscript will be omitted. This splitting of space-time allows one to 
rewrite the action (1) in the following form : 



5, 



geom 



IGttG 
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g ""{gik^f — gij^T) {Klm\n — Kin\m) 



'.^figij 



- {gijfi 



2 "-'^ 

where Kij and are respectively the extrinsic and intrinsic curvature ten- 
sors relative to the hypersurfaces mentioned above. The symbol "|" denotes 
covariant differentiation in these 3-surfaces. A^ and Ni are respectively the 
"lapse" and "shift" functions which allow one to locate the 3-surfaces with 
respect to each other. And finally, £^ denotes the Lie derivative along the 



^In the definition of the extrinsic curvature tensor, Boulware's sign conventions are 
opposite to ours 
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normal n to these hypersurfaces. 



B & L |jT6| have pointed out the existence of three main variants of the theory 
while Boulware's heuristic approach considers only two of them, i.e. when 
(3c is or not equal to zero. In this paper, we use Boulware's general formal- 
ism in the case /5c 7^ and we restrict ourselves in the next section to the 
Einstein-Hilbert Lagrangian density with a pure term. The conformally 
invariant theory (/5c = 0) will be discussed in a future paper. Using Boul- 
ware's prescriptions, the Hamiltonian form of the original action (2) can be 
written, up to surface terms, for the most general case as follows : 

Sgeom = Jd^xN {Pi, C^Q'^ + p'^ Ch9^j -nig,Q,p,P)} (3) 



with 



n{g,Q,p,P) = -2p'^P,,+acg'/^C"'''C„,k 



2 IQ-kG 
Q 
2 



R + trp2 _ p2 



+g^^P,,P+^(tr/-/) (4) 



Q'^ = (7'/2|^2ac^^o*o^ + y (5) 



2 V ^Kijj 
P = K 

where the trace and the traceless part of a tensor Aij are denoted respectively 
by A and Afj, trA^ = A'^Aij, and 



18/3c^i/2 2 ^ ' ^ 

,1/2 



iovrG 
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Q^^ is a tensor introduced in order to reduce the order of the Lagrangian as 
in Ostrogradski's method. The canonical variables are gij and Q^^ and their 
conjugate momenta are p*-' and Pij respectively. We can also write the action 
(3) in the usual Hamiltonian form, where all the constraints are manifest]^ : 



geom 



(6) 



where we have introduced the total Hamiltonian constraint Ti* as 



W 



NH + N^Hk 



(7) 



7i and Tik are respectively the super-Hamiltonian and supermomenta (or 
spatial constraints) given by (4) and 



(8) 



The constraint analysis leads to the definition of Dirac's Hamiltonian density 



(9) 



dC 



where A^'^ and act as Lagrange multipliers and H'^ = are primary 
constraints. The complete set of constraints is {'H^j,-, H^} and it is easy to see 
that are secondary constraints, i.e. that they result from the conservation 
of the primary ones, H'^. Moreover all these constraints are first-class and 
the number of physical degrees of freedom is eight in the general case. 
Dirac-Hamilton's equations are : 



9ij 



p 



dHp 

dgij 



P 



dHp 

dHp 



N„ 



H^ 



an 



D 



dHp 



A„ 



-n^" ^ 



For particular quadratic Lagrangian densities, including the cases C = R + 
PcR^ and C = dcC^^^^Cap^u, there occur new constraints. Therefore Hamil- 
tonian formalisms for distinct variants of the general theory do not use the 



^Dots denote as usual time derivatives 
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same canonical variables and are leading to distinct algebra of constraints 
1^. For our purpose, at the classical level, it is not necessary to follow B & 



L's analysis. It is sufficient to start from the general form of Boulware's ac- 
tion (3) and to impose directly the constraints that arise from the particular 
form of the Lagrangian density, when the quadratic Weyl term is neglected. 



3 Bianchi cosmological models 

With the results of the preceding section, we are now able to consider some 
cosmological applications. Our purpose is to compute the explicit expres- 
sions of the constraints and the canonical equations for Bianchi cosmological 
models in a £ = i? + ficB? theory, then to look for analytic solutions. This 
programme can be achieved with the help of REDUCE 3.5. 



The 3-metric of any diagonal Bianchi-type model can be written as follows : 



(10) 



where /i, /3_ are functions of time only and the set {uj^} is a 1-form basis 
with structure coefficients Cj;. : 

dCj' = -CL Cj^ a 
2 ^ 

In order to have a clear distinction between pure i?-squared and C-squared 
variants of the theory, it is more convenient to perform a canonical trans- 
formation from the original set of canonical variables {g^ Q, p, P} to the 
following one : 

{fi, p+, p.; n^, n+, n_; g+, q_, p+, p_, p4 
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This transformation is defined by the following relations : 

n\ = ^(2n^ + n+ + v^n_) 



g^^ = e2Mg2(/3++V3/3„) 



922 



933 



g2Atg-4/3+ 



n33 = i(n^-nH 



P, 

Q 
P 



p 



PnQr 



Ti 
3 

Ti 



V3 V3 V3 

p 

+ ^9r3 

diag (g+ + v^g_, g+ - V3g_, -2Q+] 

diag (P+ + v^P_, P+ - V3P_, -2P+) 



In terms of these new variables, the original action can be written as 

Sgeom = I {u^fi + u+(3+ + uj. + P+Q+ + p_g_ + PA - n*} 

with the 4- volume element d^Q = dt A uj^ A uj"^ A . This expression can 
be used for any diagonal Bianchi model. The explicit form of the total 
Hamiltonian constraint Ti* depends on the Bianchi model considered. In the 
theory C = R + (3cR^i we have to impose strongly the following constraints : 
Q± ~ 0. Therefore the super-Hamiltonian reduces to the following expression 



n 



R 



+ 



SttG 



e^\/*(/3+, - -e'^{P'_ + P^ - 2p2) 



(11) 



where is directly related to the usual potential in General Rela- 

tivity for the Bianchi model considered [O. 
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3.1 Bianchi IX model - pure theory 

The mixmaster cosmological model in C = theory of gravity has been 



investigated by several authors |^ |^ with the aim of studying its asymp- 
totic behaviour. 

For this model, there are no supermomenta and the super-Hamiltonian con- 
straint (11) depends on the following potential term : 

V*{P+, = ^e^^+ [cosh (4y3/?_)-l] + ^e-^^+-e-2^+ cosh (2y3/3_) (12) 

We can write Hamilton's equations resulting from (11) as follows : 

Qn 

Q± 



V3 " 



--P± 



N 

[n^ + 4P„Q„] + -^Pr^.e''^ 



V3 

N 

73 



V2 



+ P±Qn 



AtiG' 
N 



8nG 



P+e^^ 



n„ 



N 



N 



SnG 
N 

— 7='^ 



ef'V*{(3+, /?_) - 3/2e-^^(Pi + Pi - 2P, 



-2n 



Qn 



N 



V3 " 



dV* 



N dV* 
:e'^ 



Q 



We consider now the pure theory because the resulting system of canonical 
equations provides a good approximation (for large R) to the more compli- 
cated system generated by the more general quadratic Lagrangian density 
C = R + (3cR^- Due to the conservation of constraints Q± ~ 0, the canonical 
momenta P± are also constrained : 

1 

"72^ 



P±Qr. 



(13) 
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and there is no other constraint for this theory. Therefore, if Qn 7^ 0, the 
Hamiltonian constraint (11) reduces to 



-3/i 



-2/i 



from which we get the canonical equations : 



g^+-^g„y*(/3+, p_) (i4) 



ri^ = 
n± = 



AT n± 

AT 



^e-^'^Q^^-^e-^'^g„V*(/3+, /3_) 



AT 



AT 



By imposing the constraint Hr ~ and manipulating these equations, it is 
easy to get the following relation : 



(15) 



where A; is a constant of integration. When the scalar curvature is not con- 
stant, we can fix the temporal gauge by choosing ^^^P = —t which gives 
directly the following expressions for Qn{i) and N{t) : 



Qn — 2 /^cC 



N 



2k 



(16) 
(17) 



By some algebraic manipulations, it is possible to reduce the set of canonical 
equations to a differential system of three coupled equations for the physical 
variables ^{t), P±{t) : 



(18) 
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ay* 

^,, + ^, + _|e-g^t = (19) 

By direct inspection of these equations, we see that the canonical variables /i 
and (3± become uncoupled when the potential term vanishes, i.e. for Bianchi- 
type I model. 

ADM Hamiltonian 

Taking into account relations (15), (16), (17) and defining the new momen- 
tum n* , canonically conjugated to /i, by the following expression : 

n; = A; + 2P„g„ (20) 

Boulware's action 

Jd^n {U^fi + Uj+ + Uj_ + P^Q^ - NHr) 
with T-Lr given by (14), can be written as follows : 

S = Jd'Q {U^fi + n+/3+ + n_/3_ - Hadm} (21) 
where the ADM Hamiltonian density is given by the following expression : 



Ti-ADM — 777 

4kt 



(22) 



The resulting canonical equations are then written as follows 



ii = {Ul-k)l{2kt) 



P± = -U±/{2kt) 
and are equivalent to equations (18) and (19). 
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3.2 Bianchi I model - pure theory 

For Bianchi I model, there are no supermomenta and the potential term V* 
vanishes in the super-Hamiltonian constraint (14) : 

From this expression we get the following canonical equations : 



N 

N n± 



± 



2VSQr. 



n = — 



Pn 



N 
4^3 



— Qr 



and Il± = . By imposing the constraint Hr ~ and manipulating these 
equations, we get the relation (15) and we can solve the problem if we choose 
an appropriate temporal gauge. 

3.2.1 Case (^^i? not constant 

When the scalar curvature is not constant, we can fix the temporal gauge 
by choosing ^'^^R = —t which gives directly the expressions (16) and (17) 
for Qn and respectively. As stated above, the vanishing of the potential 
term V* leads to uncoupled equations (see (18) and (19)). It is now possible 
to integrate equations (19) and we get the explicit time dependence of the 
anisotropic scale functions : 

/9±(t) = -^lnt + c± (24) 

with arbitrary constants of integration c±. Equation (18) is also simplified 
and it appears to be the same equation as in Buchdahl's paper p2[. 
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It is also possible to get from the set of canonical equations the following 
first-order differential equation for : 

tt2 +172 

On the surface defined by the super-Hamiltonian constraint, its general so- 
lution coincides with the general solution of equation (18) : 



8kn 



e'^id'' + c~H~'') (26) 



where n = ^ +,9 +1 . This solution is identical to Buchdahl's solu- 



2 I fc2 



tion which was obtained from Einstein-Hilbert field equations without 

any reference to an Hamiltonian framework. 

Therefore the general form of the metric can be written : 

where we have used (17) and (24), and where /x(t) is given by (26). The 
parameters z/j are defined by 

= -^(n+ + V3n_) 



1^2 = 

2k 



l-iu^ - Vsn. 



and satisfy = ^ and J2i = 7- {4:11^ — 9). 



ADM Hamiltonian 



(n; - kf -nl- + ^p^e^^^^ 



The ADM Hamiltonian density (22) reduces to the following expression : 

(28) 

The corresponding canonical equations are then written as follows 

fi = (n;-A;)/(2H) 
$± = -n±/(2A:t) 

U± = 
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and are equivalent to equations (18) and (19) with a vanishing potential In 
particular, Buchdahl's equation mentioned above is a direct consequence of 
these equations and takes the following form : 

^' + ^ + ^^'''' = ° (29) 
It is also possible to consider another gauge than ^^•'i? = —t in order to sim- 
plify this canonical system. The new choice is ^^^R = pe*, with p = ±1, which 
can be also used for Bianchi IX model and has the advantage of including 
both classes of solutions with positive or negative scalar curvature. These 
solutions have been studied in [^. By this choice and by an appropriate 
canonical transformation defined by (20) and p* = 2p + 1, we are led to the 
following result : 

s = Jd'n {n;p* + uj^ + uj. - n^^,,} (3o) 

where the ADM Hamiltonian density is given by the following expression, 
which contains no explicit time dependence : 

1 



^ADM 



The resulting canonical equations are then written as follows 

p* = 2ni/k 



(31) 



P± = -n±/(2A;) 

1 

n± = 



The corresponding equation for p* writes as follows : 

p* - ^/e^'^* = (32) 

and its solutions exhibit the asymptotic behaviour of an evolution towards 
a de Sitter geometry when p > or a full curvature singularity when p < 



21]. The last one corresponds to the asymptotic behaviour of Buchdahl's 



solution. 
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3.2.2 Case = -4A = constant 

The variable Qn writes now Qn = —2\/?)[3c\e^^ and we have the freedom of 
choice for the function N{t). For simphcity we choose N{t) = 1. From the 
set of canonical equations and from the above expression for Qn, we get the 
following ODE : 

Qn + SXQn = (33) 

We take A > and define k = VSX. Then with a suitable definition of the 
origin of time coordinate, we get the solution for ii{t) 

g3M(t) ^ ^ 

with an arbitrary constant A. Equations (19) for P±(t) can now be solved 
easily to give 

(3±{t) = B± ln(tanA;t/2) 

with B± = (2\/3kA)^^Il±. Therefore the general form of the metric can be 
written : 

ds^ = -dt^ + sin^/^ kt Y.ta.n^''\kt/2) {dx'f (34) 

i 

where the parameters Vi are defined by 

ui = B+ + V3B^ 
1/2 = B+- V3B_ 
1/3 = -2B+ 

and satisfy = ^ and J2i = \- Buchdahl |Q, if we set i/j = rij — | 
and r = |tan(fct/2) in this metric, we get : 

ds^ = -df + cos^/^ kt J2 (dx'f (35) 

i 

where the parameters satisfy I]j = 1 and I]i = 1- If "we take the limit 
of this metric for A — > 0, it reduces to a Kasner-like metric 

ds^ = -dt^ + J2f^^ [dx'f 

i 

as expected because any vacuum solution of general relativity field equations 
is also a vacuum solution of the quadratic theory. 
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3.2.3 Isotropic CEise (Einstein-de Sitter model) 

If we look for the isotropic form of the general metric (27) we have to impose 
further constraints on the parameters Ui {ui — 1/2 — ^3). It implies that 
II-i- = and therefore Vi — and n — ^. In this case, the metric reduces 



to the following isotropic form 



ds'' = -it' + i^de + {e + ly'^dx'' + + dz'') 



(36) 



In the next section, we recover this result from the Hamiltonian formulation 
of Einstein-de Sitter model. 

4 FLRW cosmological models 



4.1 Closed FLRW model - pure theory 

We consider the isotropic metric of the closed FLRW cosmological model as 
our starting point : 



ds' = -N' de + A\t) {Cj^f + {u^f + {u'f 



with the following basis 1-forms : 



^1 ^ {i-r^f-^'^^dr 



co^^rde 



uj' — r sin 9 dcf) 



With the same prescriptions as above concerning the Hamiltonian formalism, 
we get the following action : 



where the Hamiltonian constraint is given by 



^=^-^-^^'-^^ + ^(^-3^1 (37) 
6/3c V2>A SttG \A J ^ ^ 
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and Q is defined by 



(38) 



Note that Hawking and Luttrell |P make use of another canonical variable 
Q* defined by the following expression : 

It is not difficult to perform a canonical transformation from our canonical 
variables to those used by these authors in order to write the Hamiltonian 
constraint (37) (without cosmological constant) as they do : 

n=^{Q* -Af + IlAllQ,-Q*A 

In the pure theory, the Hamiltonian constraint reduces to the following 
expression : 



V3Q 



6pc V^A 

The corresponding canonical equations are given by : 

A^Ho 



(39) 



A 

Q 
rig 



V3A 
-N 



73^2 6/3, 



' VSA 



N 



V3 + 



AQ 



4.1.1 Case ^^^R not constant 



By imposing the constraint Ti/? ~ and manipulating the canonical equa- 
tions, we get the following relation : 

AUa = QUq + k 

with an arbitrary constant k. By choosing as usual ^^^R = —t, we can reduce 
the canonical equations to the following differential equation for the scale 
function A{t) 

(40) 



AH + AA+ ^AH{12 + AH) = 



16A;2 
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As in the paper by Schmidt the integration of this equation leads to 
elhptic integrals. Note that it is not possible to avoid them by other choices 
of the lapse function. On the other end, it is easy to recover the particular 
solution A{t) = t/\^ corresponding to ^^^R = 24/t 
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4.1.2 Case = 4A = constant 

The variable Q writes now Q = 2^/3Pc^A and we have the freedom of choice 
for the function N{t). For simplicity we choose N{t) = 1. From the set 
of canonical equations and from the above expression for Q, we recover the 
exact known solution which describes a de Sitter space-time : 

A(t) = ^cosh (^^t 

If we take A = 0, the canonical equations are directly integrated and give 
the solution of Einstein's theory with incoherent radiation as source, which 
is already mentioned in Schmidt 



A{t) = 

with an arbitrary constant C. 

4.2 Einstein-de Sitter model 

The metric of the Einstein-de Sitter model writes as follows 

The Hamiltonian form of the action is 

Sgeom = I d^X [HaA + YIqQ - NU] 

where Q is already given by (38) and the super-Hamiltonian constraint writes 
as follows : 

6/3c ^flA %TiG\A J ^ ^ 

21 



In the pure R theory, it reduces to the following expression : 

The corresponding canonical equations are given by : 

NUo 



A 
Q 



V^A 



Ua = -N 



NQA 



By imposing the super-Hamiltonian constraint and manipulating the canon- 
ical equations, we find a relation among A{t) and Q{t) : 

A^ = + k2 (43) 

where k\ and /c2 are arbitrary constants of integration. We get three distinct 
cases : 

(i) If ^^^R is not constant, we choose the temporal gauge as usual by imposing 
^'^^R = —t which is equivalent to Q = —^l3cAt. Therefore equation (43) 
gives the following expression for A{t) : 



This last expression together with the canonical equations gives the function 
N{t): 



3/9c , 3/3, 



With an appropriate change of scale, expressions (44) and (45) allow us to 
write the general metric in the form (36). 

(ii) Let us note that if we set /c2 = in (43) and N{t) = 1 we get a de Sitter- 

1 /3 

like model (or anti-de Sitter) : A{t) = e^* where A = ± (y|^) • This par- 
ticular solution is recovered when we analyze the case = -4A = constant. 
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This de Sitter solution is a typical example of the relevance of the i?^-theory 
to the inflationary scenario. As a matter of fact, in the field equations, the 
terms specifically coming from the quadratic part of the Lagrangian density 
can play the role of a cosmological constant. Moreover it is well known that 
the addition of the term in the gravitational action introduces a new 
spin-0 scalar field which may act as a natural inflaton in the early universe 



(iii) When ^'^^R = 0, canonical equations are easily integrated and give, with 
an appropriate choice of gauge, the following isotropic solution : 

ds"^ = -dt^ + t {dx^ + dy"^ + dz^) 

which was found independently from a search of Groebner bases of the system 
of algebraic equations associated with power-law type solutions of the field 
equations PS . 



4.3 ADM Hamiltonian 

For the isotropic closed model considered above, it is also possible to perform 
a canonical transformation which allows one to define an ADM Hamiltonian 
density. The procedure is the same as in the anisotropic cases (subsections 
3.1 and 3.2). We set A = e°. The action can be written as follows : 

S = j d^n {li^a-TiADM] (46) 

where the ADM Hamiltonian density is given by the following expression : 

(47) 



'Hadm — -rr: 
Akt 



(Note that, for the Einstein-de Sitter model, the term proportional to e^" 
vanishes). The resulting canonical equations are then written as follows 

a = {U^-k)/{2kt) 



fl = 

8kt 
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and are equivalent to the following second-order equation 

at + a + ^ (^t'e^" + 4te^") = (48) 
which particularizes equation (18) to the isotropic regime. 

5 Conclusions 

The Hamiltonian formalism for quadratic theories of gravity, developed first 
by Boulware, has been applied to the study of the classical behaviour of 
spatially homogeneous cosmological models. We have considered here a La- 
grangian density C = R + [3cR^. 

As it is well known within the context of General Relativity the Hamilto- 
nian formalism is only consistent for the Bianchi class A models, with the 
exception of one particular Bianchi class B model, namely the type V P^]. 
In quadratic theories of gravity the situation is more complex. We have 
checked explicitly that for pure theories only Bianchi class A models can 
be dealt with the Hamiltonian formalism as described in our paper, i.e. with- 
out explicitly taking into account surface terms. However, in the case of a 
pure quadratic theory described by a Lagrangian density containing only the 
Weyl term C^^^^Cais^u, is it possible to neglect surface terms for diagonal 
metrics of class A only, with the exception of the VIq case. 
For Bianchi-type I and IX models, we have given the explicit forms of the 
super-Hamiltonian constraint, of the ADM Hamiltonian density and of the 
corresponding canonical equations. These equations are first-order and their 
compact form is well suited for analytical as well as numerical calculations. 
Therefore, their study provides an interesting alternative to the analysis of 
the usual fourth-order field equations. In the case of a pure quadratic theory 
C = l3cR^i we have solved them analytically for Bianchi I model. The solu- 
tion found is identical to Buchdhal's solution. For Bianchi-type IX model, 
we have reduced the first-order equations of the Hamiltonian system to three 
coupled second-order equations for the true physical degrees of freedom. This 
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discussion has been extended to isotropic closed as well as Einstein-de Sitter- 
type FLRW models and the explicit resolution of the corresponding canonical 
equations has enabled us to recover very easily known exact solutions. 
We also have, for Bianchi and FLRW models, performed the complete ADM 
reduction programme, in the line of Misner's pioneering work [^. The ex- 
pressions of the super-Hamiltonian constraint and of the ADM Hamiltonian 
density given here provide an adequate starting point for the quantization of 
these models, respectively in the framework of Hartle and Hawking's method 



2^ and of Misner's ADM-type quantum cosmology. 

We will tackle this problem as well as the Hamiltonian formulation of the 
conformally invariant theory (£ = acC°''^^'^Ca/3^u) in a future paper, with the 
hope of generalizing Hawking-Luttrell's work to Bianchi models. 
Another interesting possibility of dealing with the pure theory is to take 
advantage of the conformal equivalence of this theory to General Relativity 
minimally coupled to a self-interacting scalar field p9|, as it was performed 
in p for the case. However, for the conformally invariant theory, the 
situation is much more involved since it requires the introduction of a new 
spin-2 symmetric tensor besides the metric tensorial field. Accordingly the 
explicit use of this equivalence does not seem to lead, in this case, to any 
substantial simplification of the formalism in comparison to the Hamiltonian 
approach. 
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